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Abstract. We study how the singular behaviour of classical systems at bifurcations is reflected by their 
quantum counterpart. The semiclassical contributions of individual periodic orbits to trace formulae of 
Gutzwiller type are known to diverge when orbits bifurcate, a situation characteristic for systems with a 
mixed phase space. This singular behaviour is reflected by the quantum system in the semiclassical limit, in 
which the individual contributions remain valid closer to a bifurcation, while the true collective amplitude 
at the bifurcation increases with some inverse power of Planck's constant (with an exponent depending 
on the type of bifurcation). We illustrate the interplay between the two competing limits (closeness to a 
bifurcation and smallness of h) numerically for a generic dynamical system, the kicked top. 
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1 Introduction 

The semiclassical approach allows to obtain spectral infor- 
mation about quantum systems from properties of clas- 
sical periodic orbits. The most famous example of this 
quantum-classical correspondence is Gutzwiller's trace for- 
mula for completely chaotic (hyperbolic) autonomous sys- 
tems which expresses the density of states as a sum 
of contributions from periodic orbits. In such systems, the 
Lyapunov exponents A of all periodic orbits are positive, 
and their semiclassical amplitudes A cx (sinhA/2)^^ are 
finite. This is not true for systems with a mixed phase 
space, which accommodate also elliptic orbits with ampli- 
tude A oc (sina;/2)~^. The quantity uj = rwo, known as 
the stability angle, increases linearly under r repetitions 
of the orbit, and either by a suitable choice of r or of an 
external control parameter the amplitude A can become 
arbitrarily large. The contribution of an individual orbit 
eventually diverges when ll!/2tt is an integer. Normal-form 
theory ||] shows that this is precisely the condition for a 
bifurcation, the coalescence of two or more periodic orbits. 
Catastrophe theory Q further reveals that the divergence 
comes from an illegitimate stationary-phase approxima- 
tion, and provides us with uniform approximations (collec- 
tive contributions of the bifurcating orbits) that regularise 
the sing ular behaviour The ensuing true 

semiclassical amplitude is finite at ?i ^ 0, even directly 
at the bifurcation. It is clear, however, that an individ- 
ual contribution of Gutzwiller type must be valid in the 
strict semiclassical limit, as soon as one does not sit pre- 
cisely on a bifurcation (fixed distance e to a bifurcation 
in parameter space, h ^ 0). Analyticity at given U and 
varying e entails then that the true amplitude at the bi- 



furcation (e — 0) must diverge as 7i — *■ 0, with a power 
law ^ h"" as it happens to be 0. Some consequences of 
this peculiar singularity have been studied recently in the 
context of spectral fluctuations jl^. A complete solution 
of this intricate problem is emerging |l^], but it is com- 
plicated because one has to consider also more complex 
bifurcations of higher codimension which are classically 
non-generic, but are nevertheless relevant in the quantum 
realm. 

In this work we investigate the interplay of the two 
competing limits ?i — s- 0, e ^ 0, focusing on periodically 
driven systems with one degree of freedom. All results are 
relevant for autonomous systems with two degrees of free- 
dom as well. For a representative system, the kicked top, 
we find that bifurcations often interfere even when they 
are separated in phase space. This adds additional com- 
plexity to the problem at hand. We use a filtering tech- 
nique to extract contributions of bifurcating orbits and 
find that their amplitude corresponds well to the theoret- 
ical predictions. 

The paper is organised as follows: In Section ^ we de- 
scribe how the exponents for the most commonly en- 
countered bifurcations are derived. In Section ^ we present 
numerical results for the kicked top. Section ^ contains our 
conclusions. 



2 Semiclassical contributions at bifurcations 

Periodically driven systems are stroboscopically described 
by a unitary Floquet operator F. Spectral information 
about this operator is most conveniently extracted from 
the traces tr F", where n plays the role of discretised time. 
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The analogue of Gutzwiller's trace formula has been de- 
rived by Tabor pj, who found the relation 



A(x J" with 



Uq 



,- exp [US - ifi— 



(1) 



between the traces and the periodic orbits of the corre- 
sponding chaotic classical map. For convenience we denote 
here the inverse Planck's constant by = J. The sum 
is made up of all orbits of primitive (first return) period 
no with n = ngr, r an integer. The rth return of an or- 
bit is characterised by the action S = rSo, trace of the 
monodromy matrix (linearised map) M — Mq, and the 
Maslov index /i = rfiQ (which for elliptic orbits satisfies a 
slightly more involved composition law under repetitions). 

Eq. (P is valid for completely chaotic systems. For 
hyperbolic orbits, the eigenvalues of the monodromy ma- 
trix Mq are e^'*'", while elliptic orbits have unimodular 
eigenvalues e^"^". In both cases the expressions for the 
semiclassical amplitudes A oc |2 — trM| ' given in the 
introduction follow immediately. As advertised, the semi- 
classical amplitude A of an individual orbit diverges when 
Wo = 27rn/m, with n, m integers (taken relatively prime). 
The type of bifurcation depends on m, with m = 1 the 
tangent bifurcation, m = 2 the period-doubling bifurca- 
tion, m = 3 the period-tripling bifurcation, and so forth. 
Close to a bifurcation one should replace contributions of 
individual orbits by collective contributions in the trace 
formula, of the form 



A^ 



J_ 

2^ 



dl 



d0 «'(/,(/.) exp [iJ<l>(/,0)], 



(2) 



with the 'amplitude function' ^ and the 'phase function' 
(p both depending on the type of bifurcation under con- 
sideration. Here we have used canonical polar coordinates 
/, 0, which parametrise the phase space of the classical 
map as 

p—V2Js'm(l), q = V2I cos(t), (3) 

giving for the differentials dp dq ~ dl d(j). The phase func- 
tion is a local approximation to the generating functions 
S{q',p) of the classical map {q,p) {q',p'). Right at the 
bifurcation the amplitude function reduces to'P — 1, while 
the phase function is given by simple normal forms. For 
generic bifurcations we have 



(p = So ~ eq - aq^ - hp^ 
= 5*0 - - aq'^ - bp^ 
<? = 5*0 -e/- a/3/2 cos 3</> 
^ = Sq - el - aP - bp cos 40 



(to = 1), 
(to -2), 
(to = 3), 

(to 4). 



(4) 



Here e is the bifurcation parameter (bifurcations take place 
at e = 0), while Sq, a and b can be regarded as constants. 
At the bifurcation (e — 0) we can rescale the integration 
variables q, p for to = 1, 2 or / for to > 3, such that the 
combination J<P appearing in the exponent of Eq. (H) be- 
comes independent of J. What remains is a J-dependent 
prefactor in front of a J-independent integral. We find 



v = 1/6 (to = 1), 
v = 1/3 {m ~ 3), 



1^=1/4 (m = 2), 
J/ =1/2 (m > 4). 



(5) 



For £ 7^ the integral remains J-independent if one also 
rescales the bifurcation parameter according to e' = eJ^ 
[fl, with 



^ = 2/3 (to = 1), 
p = 1/3 (to = 3), 



/i=l/2 (to = 2), 
/i=l/2 (to > 4). 



(6) 



These exponents determine the semiclassical range of the 
bifurcations in parameter space. 

The case to = 3 is special in the sense that period- 
tripling bifurcations are usually accompanied by a tangent 
bifurcation, so close in parameter space that the semiclas- 
sical contribution given above looses validity for accessible 
values of J. This gives the unique opportunity to test also 
predictions for a bifurcation of higher codimension. The 
normal form is IH,nOi 



<P{I, 4>') = Sa-eI- a/3/2 3^ _ 



(7) 



The tangent bifurcation takes place at e = -^l^, while the 
period-tripling bifurcation occurs at e = 0. For e = a = 
one has to consider an integral of the form 



dl 



27T 



exp[iJ/2] oc ji/2, 



(8) 



and obtains ly = 1/2. For e, a 7^ we obtain the two 
scaling parameters /x^ = 1/2 and fia — 1/4 that charac- 
terise the semiclassical range of the bifurcation in param- 
eter space. 



3 Numerical results 

We now wish to investigate how the semiclassical singular- 
ities at bifurcations emerge for a representative dynamical 
system, the periodically kicked top |l^ , which has proven 
useful in testing semiclassical results before. The dynam- 
ics consists of a sequence of rotations and torsions, with 
Floquet operator 

F = exp ( — i ^ . I _^ j2 — iaiJzj exp (^—i/dJy^ 



X exp — 1 



2i + i 



2j + l 



J — \a2Jx 



(9) 



The angular momentum operators Jx.y,z obey the com- 
mutator relation [ Jj, Jj] — ieijkJk- Since the square of the 
angular momentum 3^ = j(j + 1) is conserved the phase 
space is the unit sphere. The role of the inverse Planck's 
constant is played hy J = j + ^, which is equal to one half 
of the Hilbert space dimension. The semiclassical limit is 
reached by sending J — + 00. We fix the rotation param- 
eters ai = 0.8, P — 1, q;2 = 0.3, and use the torsion 
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Fig. 1. The action spectrum |T*"'(5')p for the kicked top close to different types of bifurcations. For the inverse Fourier 
transformation (^) we restrict the S- integration to the intervals of width AS around the centres of the peaks Sbu, eliminating 
in this way the contributions from other periodic orbits, a) n = 1, fc = 2.6. The large peak at Sbu comes from orbits which are 
involved in a tangent bifurcation at fc ~ 2.5. b) n = 2, fc = 2.3, Qi — 1.39, close to the period-doubling bifurcation at ~ 2.1. 
c) n = 3, fc = 2, close to two period-tripling bifurcations at fc ~ 1.85 (left peak) and k ~ 1.97 (right peak), d) n — 4, k — 1.2. 



The large peak arises from orbits involved in a period-quadrupling bifurcation at k 
also involved in a period-quadrupling bifurcation, at fc ~ 1.2.) 



1.0. (The orbits of the smaller peak are 



strengths ki = k and fc2 = k/10 to control the degree of 
chaos of the classical map. The system is integrable for 
fc = and displays well-developed chaos from fc ~ 5. 

The quantum- mechanical evaluation of F is described 
in Ref. iQ . We computed the traces of the Floquet oper- 
ator and separated the contributions of different (clusters 
of) orbits by evaluating the action spectrum (a Fourier 
transformation of the trace with respect to the inverse 
Planck's constant) [la] , 



T(")(S') = 



1 



Jmax 



Jmax 



' Jmin ~1~ 1 . 



(10) 

where the difference jmax — Jmin determines the resolution 
in 5* (jmin = 1, jmax = 100). The results for parameters 
close to different types of bifurcations are shown in Fig. |l]. 
The contribution at given j of orbits pertaining to a given 



peak can be obtained by an inverse Fourier transforma- 
tion, 



Aoc J dsl ^ trF"(j')e-'(^'+^'^ e'(^+^)^ 
= 2 X; trF"(/)e'(^-^')^B. smjU -f)AS/2 



i'=jmin 



(11) 



where the integral over actions S is restricted to an inter- 
val AS around the centre S'sif of the peak. This eliminates 
contributions of other periodic orbits. 

It is convenient to tune the control parameter fc slightly 
away from the bifurcation to the value that maximises 
the contribution of the bifurcating orbits to trF". The 
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Fig. 2. Logarithmic plots versus inverse Planck's constant J 
of the maximal (in parameter space) contributions | A| of orbits 
involved in the bifurcations of Fig. |l[ calculated by Eq. (|ll|. a) 
Tangent bifurcation. The average (dotted line) gives the expo- 
nent 1/ ~ 0.1866. b) Period-doubling bifurcation, f ~ 0.2636. c) 
Period-tripling bifurcation of higher codimension, i/ ~ 0.5327. 
d) Period-quadrupling bifurcation, v ~ 0.5734. 



parameter k of the maximum approaches the true bifur- 
cation point with the exponent fj,, Eq. (^. The maximal 
contribution is of the same order of magnitude as the 
contribution at the bifurcation, but it is less sensitive to 
changes in the parameters. Most importantly, this proce- 
dure does not require any classical information (like the 
precise parameter value of the bifurcation) and is hence 
genuinely quantum mechanical. We extract the exponents 
p from logarithmic plots of the maximal \A\ versus J, 
shown in Fig. ^. In all cases we find good agreement with 
the theoretical predictions. For a tangent bifurcation at 
k ~ 2.5 the observed exponent is ~ 0.1866 (theoreti- 
cally, u = 1/6). Two different period-doubling bifurcations 
appear at k ~ 2.8 and produce overlapping peaks in the 
action spectrum. We separated them by changing ai to 
ai = 1.39, moving in that way one of the period-doubling 
bifurcations to k ~ 2.1. The exponent for this bifurcation 
is ~ 0.2636 (theoretically, v = 1/4). Back to the original 
value ai — 0.8, we find for the period-quadrupling bifur- 
cation at fc ~ 1.0 the exponent i' ~ 0.5734 (theoretically, 
z.= l/2). 

As mentioned above, for all known dynamical systems 
period-tripling bifurcations are typically accompanied by 
a tangent bifurcation, so close in parameter space that 
one has to treat the situation as a bifurcation of higher 
codimension. For the kicked top, an angular momentum 
of about J ~ 10^ would be needed for separating the orbits 
in the 'period-tripling-|-tangent' bifurcation at fc ~ 1.85. 
The same is true for a similar sequence of bifurcations at 
k ~ 1.97. For the much smaller values of J that we use 
here we hence have the unique opportunity to test the 
exponent for a case of higher codimension. As before, the 
result 1/ ~ 0.5327 (for the bifurcations at fc ~ 1.85) is close 
to the theoretical expectation v = 1/2. 



4 Conclusions 

We have studied the asymptotic behaviour for ?i of 
periodic-orbit contributions to semiclassical trace formu- 
lae, around points in parameter space where orbits bi- 
furcate. For the most common types of bifurcations the 
theoretically predicted power-law divergence oc h^'^ was 
tested numerically for a representative dynamical system, 
the kicked top, giving good agreement for the exponents 
i>. 

In the semiclassical limit the contribution of non-bi- 
furcating orbits reaches a constant value \A\ = 0{?i^), 
corresponding to i/ = 0. It follows from Eq. (||) that 
the exponent for bifurcating orbits falls into the range 
< J' < 1. As a consequence, the semiclassical contribu- 
tion of bifurcating orbits is dominant when parameters are 
close enough to the bifurcation point. On first sight this 
seems to require a careful tuning of the parameters. From 
the perspective of spectral statistics, however, a careful 
tuning often turns out to be unnecessary ||ll|,|l2j: The pe- 
riod n of orbits that contribute to the spectrum increases 
in the semiclassical limit as well, and one enters a com- 
petition between the weight of bifurcations in parameter 
space (given by the exponents v and and the prolifer- 
ation of their number with increasing n. Some quantities 
are dominated by bifurcating orbits even when the pro- 
liferation is not taken into account. The final outcome of 
this competition is not clear at the moment and certainly 
deserves further investigation. 
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